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I review the current status of lattice calculations of heavy quark quantities. Particular emphasis is placed on 
leptonic and semileptonic decay matrix elements. 



1. INTRODUCTION 

B and D meson hadronic matrix elements play 
an important role in determinations of the CKM 
matrix elements and overconstraining the unitar- 
ity triangle of the Standard Model. Form fac- 
tors, decay constants and bag parameters are de- 
termined from these hadronic matrix elements 
which can be calculated directly on the lattice. 
B factories (BaBar and Belle) and Charm fac- 
tories (CLEO-c) will reduce experimental uncer- 
tainties in measurements of decay process, leav- 
ing the theoretical calculations as the dominant 
uncertainty. Although the focus of lattice heavy 
quark calculations has traditionally been on B 
physics, new experiments such as CLEO-c will 
determine fo^^f and semileptonic branching frac- 
tions to a few percent, testing lattice results for 
these quantities. 

In this talk I will report on progress made in the 
last year in calculations of heavy quark quanti- 
ties. Many quantities have now been determined 
by different groups using different heavy quark 
methods. The talk is organised as follows. I be- 
gin by discussing the status of quenched and un- 
quenched calculations of the leptonic decay con- 
stants. Progress in determinations of the B pa- 
rameters, Bb^j^^^, which together with are 
relevant to determinations of \Vtd\ and \Vts\ and 
to Bg width and b hadron lifetime measurements, 
is then reviewed. Turning to semileptonic decays, 
I will focus on B ^ -kIv and a new method to de- 
termine the zero-recoil form factor of B D*lv. 
A new calculation of the h quark mass and an 
outstanding puzzle in lattice calculations of spin 
splittings are briefly discussed. 



1.1. Heavy quark methods 

The calculations described here employ three 
different heavy quark methods, each of which 
has different systematic errors. The UKQCD 
and APE groups use a nonperturbatively 0{a)- 
improved Sheikholeslami-Wohlert (SW) action 
(NPSW) at quark masses about that of charm. 
Results are extrapolated to the h quark guided 
by heavy quark effective theory (HQET). A large 
source of uncertainty comes from this extrapola- 
tion. It is well known that the nonperturbative 
C'(a)-improvement scheme breaks down in the 
limit mg oo. This was discussed last year by 
Bernard and an explicit example was given at 
this conference by Kurth and Sommer They 
have shown that the continuum limit (a 0) 
must be taken before an extrapolation in heavy 
mass is attempted. 

In nonrelativistic QCD (NRQCD) g the heavy 
quark is assumed to be nonrelativistic ((w/c)^ ~ 
(0.3GeV/5.0GeV)^). Relativistic momenta are 
excluded by introducing a finite cut-off such that 
p sa mqv rriQ. Then, IpI/uiq ^ 1 and the 
QCD Lagrangian can be expanded in powers of 
l/iriQ. This has been successful for B physics, 
where the quark mass is large, so the expansion 
is convergent. The theory is nonrenormalisable 
but at finite lattice spacing, lattice artefacts can 
be removed by including higher orders in l/mg 
and a. 

The Fermilab approach (FNAL) |^, as used 
in current heavy quark calculations, is a re- 
interpretation of the SW action which identifies 
and correctly renormalises nonrelativistic opera- 
tors present in the SW action. Discretisation er- 
rors are then 0{aAQCD) and not 0{amQ). The 
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existence of a continuum limit means a contin- 
uum extrapolation is possible. An implementa- 
tion of the full approach, for C'(o)-improvement, 
has been done by the Fermilab group Q . 

2. LEPTONIC DECAY CONSTANTS 

The B meson decay constant is an input in 
determinations of the CKM matrix element \Vtd\ 
through — B^ mixing. The mixing, Arud, is 
known quite precisely and the dominant uncer- 
tainty comes from Jb^JBe^- 

2.1. The quenched approximation 

A compilation of recent determinations of Jb 
and f]j is shown in Figs |l| and ^ The calcula- 
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Figure 2. Recent determinations of fo in the 
quenched approximation. 
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Figure 1. Recent determinations of fs in the 
quenched approximation. * indicates results up- 
dated from Lattice'OO. 

tions are ordered by the heavy quark action used. 
The solid bands indicate my average for the par- 
ticular heavy quark action used and values of Jb 
obtained for a given heavy quark action can be 
seen to be in good agreement. In addition there 
is a clear overlap between results from the differ- 
ent heavy quark actions. 

The most recent calculations of in Fig || also 
show there is broad general agreement between 



groups and between the two treatments of heavy 
quarks. 

In 2001 the CP-PACS and UKQCD collabo- 
rations both reported updates to their quenched 
determinations of /s, in Refs [|3| and |16 1^ 
respectively. In this CP-PACS calculation an 
NRQCD action, corrected to ©(l/M), is used 
and /s is been determined in both the quenched 
approximation and using unquenched gauge con- 
figurations. The unquenched result is unchanged 
and is included in Table ^ The central value 
in the quenched approximation is also unchanged 
from that discussed at Lattice'OO by Bernard 
However, the systematic error has been revised 
upwards to reflect the observed large discretisa- 
tion effects. The estimate of the continuum value 
is now Jb = 191 ± 4 ± 27 MeV and Jb, = 
220 ± 4 ± 31 MeV, where the first error is sta- 
tistical and the second is the uncertainty due to 
discretisation. There are additional uncertainties 
in excess of 30% from setting the lattice scale and 
3% in from ambiguities in determining the 
strange quark mass. 

Although the central value is in reasonable 
agreement with results from other groups such 
large discretisation errors have not been observed 
previously. This scaling analysis is done us- 
ing fairly coarse lattices, 1.017GeV < < 
1.743GeV but, nonetheless, both GLOK "and 
JLQCD who also use an NRQCD action corrected 
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to 0{1/M) and take operator mixing fully into 
account, as in the CP-PACS calculation, see no 
such scaling violations. The GLOK group were 
the first to show that including the ©(a^a) dis- 
cretisation correction to the heavy-light axial cur- 
rent leads to a significant decrease in scaling vi- 
olations of fs- Indeed both GLOK and JLQCD 
find quite mild dependence on the lattice spacing, 
in contrast to CP-PACS. The calculations differ 
in the choice of gauge action. CP-PACS use an 
Iwasaki RGI action which although it has been 
shown to have small scaling violations, the decay 
constants seem to be an exception ■ 

It is important to understand the reason for 
the observed scaling violations, especially in a 
quenched calculation, where a motivation for the 
approximation is that it allows us to control and 
understand other systematic uncertainties. 

The systematic errors in the calculation of de- 
cay constants by Lellouch and Lin in Ref |l6| have 
been revised to include the effects of using differ- 
ent physical quantities to determine arris. The 
changes involved are very small. 

The results of the UKQCD Collaboration, re- 
ported in 1999 in Ref and discussed last 
year in Ref [Q have been substantially revised. 
The new analysis was described by Maynard at 
this conference jl^. UKQCD determine /b us- 
ing a non-perturbatively O(a)-improved action, 
extrapolated from charm to the bottom quark 
mass. The first change is in the current and mass 
improvement and renormalisation coefficients. A 
consistent set of nonperturbative coefficients, as 
determined by Bhattacharya et al [|l9| is used. 
Secondly, the scale is set by the pion decay con- 
stant rather than ro as previously. This is a better 
choice since rg is poorly determined experimen- 
tally. Finally, an extended discussion and analy- 
sis of the systematic errors and in particular the 
error in the heavy mass extrapolation is included. 
The mass-dependent normalisations proposed in 
the Fermilab formalism Q| and by Bernard [Q 
are compared to the nonperturbative normalisa- 
tion which is used to determine the central value. 
The effect of 0{a^m^) lattice artefacts at the 
charm scale on the extrapolation to the bottom 
quark mass is studied by simultaneously fitting 



data from two lattice spacings to 

Wa) = 7(1 + ^ + ^ 

+e{aMf + ^{aMf) . (1) 

This is the usual HQET scaling relation with ad- 
ditional mass-dependent terms. From the fit a 
so-called "quasi-continuum" result is determined 
and compared to fits at individual lattice spac- 
ings, the difference being taken as a measure of 
the effect of C'(a^mQ) lattice artefacts. While the 
method proposed here is a reasonable one, the size 
of the error envelope at the B meson mass from a 
fit to Eqn ^ would be a more realistic estimate of 
the uncertainty in this extrapolation. This would 
almost certainly be larger than the error currently 
quoted by UKQCD. 

I conclude this section with a description of a 
calculation by Davies et al of /b at rest and at 
non-zero momentum, using NRQCD. The moti- 
vation for such a study comes from the need to un- 
derstand momentum-dependent errors in semilep- 
tonic decays, for which fs offers an easy place to 
start. The matrix element {0\A^\B{p)) — fsPfj. 
is studied with temporal, Aq and spatial, Ak cur- 
rents. A number of different smearings are inves- 
tigated and narrow smearing of heavy quarks is 
found to be optimal for moving B mesons. Con- 
strained curve fitting |^ gives reliable results at 
momenta up to (pa)^ = 16, in units (27r/i) and 
good agreement is found for fs at zero and non- 
zero momentum. Finally, the discrepancy dis- 
cussed in Ref between fs from Aq and Ak 
is resolved by including a correct power counting 
and appropriate normalisation of the different 
in the matrix element. 

2.2. Unquenching 

The last year has seen some new calculations 
of decay constants with dynamical quarks. These 
and other recent calculations are summarised in 
Table 0. As in the quenched case these calcula- 
tions are performed with a number of heavy quark 
actions. I will discuss the results marked as "new" 
in Table |^. Details of the other calculations can 
be found in the references quoted and in Ref Q|. 

Yamada [|2^, for the JLQCD Collaboration, 
presented preliminary results from a calculation 
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Table 1 

Leptonic decay constants with Nf = 2. A * indicates a new result. 



Group 


/B(MeV) /^(MeV) 


Collins99 |24|] 
MILC'OO*!^ 
MILC'Ol* {Nf = 2 + 1) |27| 
CP-PACS'OO(FNAL) [|| 
CP-PACS'OO(NR) @ 
JLQCD* p5ll 


186(5)(25)(+g'') ~ 1.26 

191(6)(l?^)(tji) ^1.10 215(5)(±i^)(t«) ^1.08 
1.23(3)(11) 

208(10)(29) 1.11(6) 225(14)(40) 1.03(6) 
204(8) (29) (l;5^) 1.10(5) 
190(14)(7) ~ 1.14 



of /s(^) using NRQCD, including all \/M correc- 
tions. A nonperturbatively C'(a)-iniproved Wil- 
son action for the sea and valence hght quarks 
and a plaquette gauge action are used. They 
compare values of Jb in the quenched approxima- 
tion using a perturbative (tadpole-improved) and 
nonperturbative value of csw- The difference is 
found to be negligible. Using nonperturbative im- 
provement coefficients, the quenched value of Jb 
is compared to the two flavour result {Nf = 2). 
/b is larger in the unquenched theory but the 
size of the difference depends on the type of va- 
lence chiral extrapolation used. Linear, quadratic 
and log fits were compared and the latter two 
yield good x^/^df- The log fits are motivated by 
chiral perturbation theory ||2^ and to my knowl- 
edge this fit form has not been explored by other 
groups, for this particular quantity. The results 
from JLQCD indicate a careful analysis of the chi- 
ral behaviour is warranted. For their preliminary 
result, in Table ^ the central value is an average 
of the quadratic and log fits and the second er- 
ror is the difference from the linear extrapolation. 
They are gathering more statistics for a detailed 
study of the chiral extrapolation. The group also 
presented results for the B parameters which 1 
will discuss below. 

The "MILC'OO" result was discussed in some 
detail in Ref . I refer the reader there for more 
details and turn to the newest result from MILC. 
Preliminary results for fs and ratios were pre- 
sented at this conference by Bernard |27|. An 
improved Kogut-Susskind action ( "Asqtad" ) with 
Nf = 2 + 1 was used. For heavy quarks they use 
a tadpole-improved SW action with the Fermilab 
nonrelativistic interpretation. The dependence of 



/b on dynamical {rriT^/mp)'^ is investigated with 
constant and linear fits. The valence chiral ex- 
trapolation is seen to be well controlled with the 
difference between linear and quadratic fits be- 
ing quite small, in contrast to the JLQCD result. 
Since both calculations are preliminary, definite 
conclusions are premature and further investiga- 
tion is required. The currents have not been nor- 
malised and so 1 do not include the value of /s in 
Table |l|. In a ratio of decay constants the normal- 
isation factors largely cancel, as do many system- 
atic uncertainties. MILC find a mild dependence 
on both Nf and a^^ in the ratio fB^I Jb- 

Unquenched simulations (with Nf = 2) of D 
meson decay constants have been done by CP- 
PACS and MILC as reported in Table [|. Both 
groups use the Fermilab nonrelativistic interpre- 
tation for the heavy quarks and an improved Wil- 
son action for light quarks, although with differ- 
ent gauge actions. The results agree within the 
quoted uncertainties. The results at three lattice 
spacings do not show scaling and CP-PACS note 
that Jb may in fact be smaller than their reported 
value due to discretisation effects. They are able 
to conclude however that the dynamical effect for 
D mesons is smaller than for B mesons. This is 
supported by the MILC data. 

It is a difficult task to combine results given 
the very different systematic errors in different 
calculations. Nonetheless, I believe the estimates 
below represent a summary of the current status. 



Quenched 
Jb - 173(23) MeV 
Jb, = 200(20) MeV 
Jd = 203(14) MeV 
/d, = 230(14) MeV 



Unquenched 
/b = 198(30) MeV 
fB, = 230(30) MeV 
/d = 226(15) MeV 
Jd, = 250(30) MeV 
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Id, I fo 



1.15(3) 
1.12(2) 



/s, / fs ■ 



1.16(5) 
1.12(4) 



Of these, only has been measured experi- 
mentally. A recent average is Jd^ = 280 ± 48 
MeV |^8|, in agreement with the average above. 

3. B PARAMETERS 

In comparison to lattice calculations of B 
parameters have not received as much attention. 
There is however, progress to report this year. 

A new calculation of AB — operators by APE 
is described in Ref These matrix elements 

are important in lifetime ratios of B mesons. 
Reyes ^] presented details of a quenched cal- 
culation by the APE collaboration of the matrix 
elements of AB = 2 operators. A more detailed 
description has recently appeared in Rcf |3l) . Re- 
sults from a relativistic simulation (at the charm 
quark mass) are combined with results in the 
static limit and interpolated to the b quark mass. 
The operators in the relativistic and static theo- 
ries are matched at NLO, including the anoma- 
lous dimension matrix up to two loops and the 
one-loop matching for the AB = 2 operators. 
This is the first time such a consistent match- 
ing has been done and since the AB = 2 op- 
erators are determined by interpolation rather 
than extrapolation the systematic error is sig- 
nificantly reduced. The result is BB^imi,) ~ 

0.87(4) (3) (0) (it), where the errors are statis- 
tical, systematic errors in the static calculation, 
error in the renormalisation in QCD and a com- 
bined error due to the uncertainty in a^^ and in 
the improvement coefficient of the axial current. 
Not included here is a full estimate of the dis- 
cretisation error although it is discussed in the 
text. The operators are not C'(a)-improved and 
the static and SW calculations are done at differ- 
ent lattice spacings. The discretisation error due 
to the unimproved current is included in the third 
error bar, any remaining a-dependence is difhcult 
to estimate with just one lattice spacing. 

JLQCD H use NRQCD determine B parame- 
ters from the quenched and unquenched data sets, 
already described in Section 2^ . They find a neg- 
ligible difference in the two results. The matrix 



elements are renormalised at the one-loop level 
in perturbation theory. The light quark mass de- 
pendence is mild and no significant dependence 
on the choice of fit function is observed. The de- 
pendence on the quantity used to set the scale and 
the determination of the strange quark mass are 
included in the error budget. Their result, with 
iV/ = 2 is BsAmb) = 0.872(39)(19)(4). The er- 
rors are statistical, the difference from four meth- 
ods to determine Bb and the chiral extrapolation. 
Power counting leads to an additional 8% system- 
atic error. Sharpe has estimated |32j the effect of 
the quenched approximation, from quenched chi- 
ral perturbation theory, to be ~ 10%. This cal- 
culation suggests that this 10% may be an upper 
bound on the effect of unquenching but as this is 
a preliminary result I think it is best to leave the 
10% uncertainty and await the final result and 
subsequent confirmation by other groups. 

Taking these two new (albeit preliminary) re- 
sults into account I estimate the renormalisation 
group invariant B parameters to be 



Bb^ 
Bb^ I Bb^ 



1.30(12)(< 10%) 
1.34(10)(< 10%) 
1.03(3) 



fB^VB^JfBy/B^, = 1.15(6) i 



The second, asymmetric error on ^ takes into ac- 
count the preliminary result from JLQCD. This 
suggests that in the unquenched calculation of fs 
chiral logs are important in the chiral extrapola- 
tion, fs^ is uneffected so the value of ^ in future 
may be larger than presently quoted. 

4. SEMILEPTONIC DECAYS 

A determination of |Kif,| and |Vc&| can be made 
by combining experimental measurements of ex- 
clusive branching fractions of B —^ it, plv and 
B D'^*Hv with lattice calculations of the 
hadronic matrix elements (or corresponding form 
factors). In practice, lattice calculations of these 
matrix elements, with experimentally favoured 
kinematics present various problems. I will re- 
port on some recent progress in calculations of 
B(D) TT.plv and B D^*^^. Preliminary re- 
sults, from UKQCD for the decay B plv were 



6 



presented at this conference by Gill ||37|| . 

5. B -> TTliy 

The hadronic matrix elements parameterising 
the decay of heavy-light to light mesons are not 
protected by flavour symmetries and hence \Vub\ 
is dominated by theoretical uncertainties. 

\Vub\ is determined from the decay rate 



dp 



G%\Vut\' 2mBpV+iE)\' 



247r3 



E 



(2) 



where the form factors, f+{E) and fo{E) param- 
eterise the matrix elements in the standard way 



{7r{p^)\V^\BipB)) 



,2 1 M 



f+{E)[pB+Pn 



-r 



(3) 



I will focus on four recent calculations by 
UKQCD jp, APE g, FNAL || and 
JLQCD |§T^ UKQCD and APE extrapolate 
from quark masses around charm to the bottom 
quark mass. The Fermilab group use the Fermi- 
lab formalism for heavy quarks and JLQCD use 
NRQCD. Thus it is possible to investigate the 
agreement of different analyses and treatments of 
heavy quarks. Table ^ lists some details of these 
analyses. A number of improvements, compared 
to previous calculations are included. All groups 
make the chiral extrapolation and the renormal- 
isation of matrix elements is done nonperturba- 
tively or at one loop in perturbation theory. The 
lattice spacing dependence is studied by the Fer- 
milab group who also simulate directly at the b 
quark mass so there is no uncertainty due to ex- 
trapolation of heavy quark mass in the calcula- 
tion. This is, of course, also true of the JLQCD 
calculation. 

Fig H shows the differential decay rate as a func- 



tion of the recoil momentum, 



q 



'2mBE- 



TO^, from the four groups listed above. The error 
bars are the statistical and systematic errors in 
each calculation. The data show general agree- 
ment between the four groups, within the quoted 
uncertainties. The major shortcoming of all lat- 
tice calculations can also be seen from this plot: 
the bulk of the experimental data is at close to 
zero. To avoid a model-dependent extrapolation 



in q^ new techniques are needed to increase the 
range of momenta, at currently accessible lattice 
spacings. Alternatively new experiments may be 
able to extend the current kinematic range for 
this decay. 



□ APE 
O Fnal 
O JLQCD 
< UKQCD 




Figure 3. Comparison of differential decay rates. 



5.1. Fermilab and JLQCD 

Rather than use the paramctcrisation of Eqn ^ 
both Fermilab and JLQCD use forms motivated 
by HQET. 



(^(p,)|V^|i?(pB)> 



V2^[v^fl\{E) 



where v = PB/mB and p± = p-n — Ev. The first 
parameterisation is used by Fermilab and the sec- 
ond by JLQCD. The traditional form factors, /+ 
and fo are easily obtained from /|| and f± or /i 
and /2. These are useful quantities to focus on. 
Considering the , E limit yields 



/II 



f± = 



V2U 



V2U 



-gBB'TT- 



2mB 



(4) 



(5) 



through ©(I/toq) in the heavy-quark expansion. 
Similar expressions can be written down for /(i,2) ■ 
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Table 2 



Overview of the B ttIv analyses. PT is perturbation theory and NP is nonperturbative. 





APE 


UKQCD 


FNAL 


JLQCD 


simulate at nib 


no 


no 


yes 


yes 


quark mass 




M,est 


Mkinotic 


Mkinctic 


chiral extrap. 


yes 


yes 


yes 


yes 


lattice spacings 


1 


1 


3 


1 


matching 


NP at = 


NP at ruq = 


mostly NP-.y/Zv^'-Zytb 


PT 




includes KLM term 


includes KLM term 


PT: radiative corrections 
at 1-loop 


1-loop 



In addition, these form factors have simple de- 
scriptions in HQET and so are natural in the 
Fermilab and NRQCD approaches. Finally, they 
emerge directly from the lattice calculation and 
can be analysed directly, forming /(+.o) at the 
end. 

In Ref a discrepancy between one of the 
Fermilab and JLQCD form factors was discussed. 
Fig ^ shows this discrepancy in fi + f2- This form 
factor yields /o which does not contribute to the 
rate. However, one would expect the NRQCD 
and Fermilab approaches to be in agreement es- 
pecially as the simulation parameters are also 
similar. In fact the functional form of the chi- 
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Figure 4. Comparison of JLQCD and Fermilab 
form factors 

ral extrapolation and the range of light quark 
masses has a significant effect on the final val- 
ues. JLQCD use a linear extrapolation with a 
range of light quark masses 0.034 < amo < 0.072 



and Fermilab determine the central value with 
from a quadratic fit with light quarks in the range 
0.024 < amo < 0.059. With this dependence on 
the light quark masses and the functional form of 

the chiral extrapolation, the ^ ^ifj % uncertainty 
due to chiral extrapolation assigned to dV /dp in 
the Fermilab analysis seems prudent. It is in 
fact the largest systematic error in that analysis. 
Lighter quark masses and higher statistics would 
certainly help to reduce this uncertainty. 

5.2. Heavy mass dependence 

Turning now to the heavy mass dependence. 
This has been studied by JLQCD and APE al- 
though not in the same region of heavy quark 
mass. JLQCD simulate at four heavy quark 
masses around the b quark. Plotting the quanti- 
ties = (a,(Mp)/a«(Ms))-2/"/+/\/Mp and 
$0 = (a.(Mp)/a,(MB))-2/ii/oVMp at fixed E 
they see no dependence on the heavy mass. This 
analysis is restricted to the bottom quark regime. 
On the other hand, APE find a negative slope 
when plotting the same quantities versus 1/Mp 
with heavy quarks at around the charm mass. 
These differences in slope can be seen in Fig- 
ure H The Fermilab nonrelativistic approach is 
valid at arbitrary quark mass. Fig ^ compares the 
JLQCD, APE and Fermilab data where only Fer- 
milab have simulated directly at both the charm 
and bottom quark masses. The figure reveals that 
the Fermilab data agree with JLQCD at Mb but 
show slope as 1/Mp 1/Md- The slope of the 
mass dependence is gentler than that observed by 
APE agreeing more with a linear extrapolation of 
the APE data to Mb ■ 

I make some final remarks in this section about 
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Figure 5. Mass dependence of JLQCD (o), APE 
(o) and Fermilab (□) form factors. 



the soft pion theorem (SPT). The relation appHes 
in the Hmit TOtt , i? as given by Eqn ^. The 
first remark is that the hmit in which this relation 
holds implies /o should be evaluated at rn^ rather 
than g^^^ = m% - ml. 

The question of whether or not lattice calcu- 
lations satisfy this relation has been discussed in 
previous reviews 0j3^ . In the analyses discussed 
here both APE and UKQCD report that the rela- 
tion in Eqn^ holds. In the JLQCD and Fermilab 
results this relation does not hold. It is not at 
all clear what the reasons for this are. The APE 
group use the soft pion theorem as a criterion 
for choosing the functional form (quadratic rather 
than linear) of the heavy mass extrapolation so it 
is perhaps not surprising then that the relation 
holds. UKQCD also choose a quadratic form for 
the heavy mass extrapolation and there is a de- 
pendence on the pole form introduced when inter- 
polation to a common set of q^. JLQCD observe 



that a , 



■ term in the chiral extrapolation raises 



the value of the soft pion limit. Thus the problem 
is dependent on the form of the extrapolation and 
the interpolation to fixed E. In addition, Onogi 
has pointed out ||3^ that the perturbative uncer- 
tainty in the renormalisation coefficients of the 
heavy-light currents may be contributing to the 
SPT discrepancy. This year Onogi |Q reported 
on a nonperturbative determination of the ratio, 
Za/Zv with static heavy quarks and SW light 



quarks. The result brings the two sides of Eqn ^ 
into better agreement, but does not explain com- 
pletely the difference. 

Finally, the Fermilab group note that while the 
lattice data for /o do not satisfy Eqn ^ the exper- 
imental rate at this point goes to zero and the 
lattice systematics (mostly due to chiral extrap- 
olation) increase at the soft pion limit. For these 
reasons cuts are imposed on the range of energy 
considered such that E > 0.424GeV. Although 
the phenomenologically interesting region is far 
from g^Qj. the reasons for the violations of the 
SPT are important. It is a good indication of the 
control of systematic errors. 

The sources and severity of the different sys- 
tematic errors in the calculation of differential de- 
cay rates and their contribution to the theoretical 
uncertainty in \Vub\ were tabulated in the Fermi- 
lab analysis. The details are in Table |. Tb IS 

Table 3 

Sources of systematic error, in percent, in 
B{D) ttIv from the Fermilab analysis psf. 



Error 


Tb 


\Vub\ 


Td 




statistical 


+27 
-9 


+ 14 
-5 


+ 17 
-8 


+9 
-4 


excited state 


±6 


±3 


±6 


±3 


p extrapol" 


±10 


±5 


±9 


±5 


rUq extrapol" 


+ 16 
-22 


+8 
-11 


+3 
-18 


+2 

-9 


adjusting toq 


±6 


±3 


±2 


±1 


HQET matching 


±10 


±5 


±10 


±5 


a dependence 


+ 16 

-3 


+8 
-2 


+23 

-6 


+11 

-3 


definition of a 


±11 


±6 


±4 


±2 


Total systematic 


30 


15 


+28 

-24 


+ 14 
-12 


Total (stat(g)sys) 


+40 
-31 


+20 
-16 


+32 
-26 


+ 16 
-13 



defined from Eqn | as Ts = 2mBP^\f+{E)\^ / E 
and similarly for Tb- The D ^ t:Iv calculation is 
discussed in Section ^. 

6. D ^{Kyv 

This decay has traditionally received less atten- 
tion than the B nlv. With the planned charm 
factories at CLEO-c and elsewhere it is increas- 
ingly important to focus on D meson decays. This 
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is a region of quark mass where Fermilab's non- 
relativistic interpretation and the C'(a)-improved 
approach should be in agreement. A nice advan- 
tage in this decay is that lattice calculations can 
reach the entire kinematical range with no ex- 
trapolation in . This removes what is one of 
the major problems with the B meson decay. 

Both APE and Fermilab have reported results 
for the decay D ttIv. In addition the APE 
group presented their results for D — > Klv. Pre- 
liminary results for this transition were presented 
by the Fermilab group in Ref pd|. 

As already indicated by Fig |5| the results of 
the Fermilab and APE analysis do not agree at 
the D meson. The difference is even greater 
when one looks at the differential decay rate since 
(dT/dq'^) (X There are a number of 

differences in both the action and improvement 
coefficients and the analysis of the groups. The 
Fermilab group define the quark mass from the 
kinetic mass rather than the rest mass, as APE 
does. In lattice calculations of heavy quark sys- 
tems the kinetic and rest masses do not agree. 
The kinetic mass describes the dynamics of a non- 
relativistic system and this choice is advocated 
in the Fermilab formalism (as well, of course, 
as in NRQCD). However, at masses around that 
of charm the difference between rest and kinetic 
mass is not as great and the choice of one defini- 
tion rather than the other cannot, I believe, ex- 
plain the lack of agreement. Further differences 
arise in the current renormalistation. In the Fer- 
milab scheme the bulk of this is done nonper- 
turbatively in a fully mass-dependent way. The 
remaining matching from lattice HQET to contin- 
uum HQET is done at one-loop level in pertur- 
bation theory. APE use a nonperturbative deter- 
mination of the matching coefficients for massless 
quarks with the so-called KLM term to correct for 
0{ma) effects. In addition there are differences 
such as the use of a quadratic (Fermilab) and a 
linear (APE) chiral extrapolation. It may be that 
these differences combine to produce the large dis- 
crepancy seen in the form factors. With experi- 
ments soon to test lattice predictions it is impor- 
tant to understand what is going on. I would 
also like to note that the UKQCD collaboration 
do have preliminary results. Since they are not 



final I have not included them in this discussion 
although interestingly they lie between the Fer- 
milab and APE results. 

7. B ^ D^*Hv 

New results were presented in two talks at this 
conference. Simone Q| reported on work by the 
Fermilab group to measure the form factors for 
this decay at zero recoil. Lacagnina reported 
preliminary results from UKQCD for the shape 
of the Isgur-Wise function determined from B 
D*lv and B — > Dlv. I will focus the Fermilab 
results since they were (almost) final at the time 
of the conference and the calculation has recently 
been completed and appeared in Ref [Q. 

At the zero recoil point the simple relation 
Tb->d* (1) — (1) holds and heavy quark sym- 
metry constrains to have the heavy quark 
expansion 



1 



h 



21, 



h 



(2toc 



2mr2mh 



(2mb)2 



.(6) 



In Eqn |^ the radiative correction, tja is known 
to two-loop level from a calculation by Czarnecki 
and Melnikov The three Is are hadronic ma- 
trix elements of the HQET and are calculable in 
lattice QCD. 

To determine these matrix elements, for the 
first time in a lattice calculation, a new method 
is introduced. It exploits an idea developed in a 
previous paper where double ratios of matrix ele- 
ments were used to determine the form factors of 
B Dlv at zero recoil [Q. Three such double 
ratios are employed to determine ly^lp and I a- 



TIa^ 



{D\n%\B){B\b-i'^c\D) 
{D\c-i^c\D){B\}ri%\B) 
{D*\c^%\i3*){i3*\b-i'^c\D*) 
{D*\c-i^c\D*){B*\}ri'^b\B*) 
{D*\c-fj-f^b\B){B*\lri^-i^c\D 



{D*\c^,^5c\D){B*\b^,-^5b\B 



-\hM{l)? 



where |A^,(l)p is defined by ^^^^r^ili^^C^l;} ■ 

The matrix elements lp,lv and I a are obtained 
from the heavy quark expansions of /i+(l), /ii(l) 
and /lAi(l) respectively. The required form 
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factor, /i^j^(l) can then be reconstructed using 
Eqn ||. The Fermilab group report a final result 



Tb^D'UI) = 0.9130 



+0.0238+0.0171 
-0.0173-0.0302 



(7) 



where the first error is statistical and fitting and 
the second systematic, resulting from matching 
lattice gauge theory and HQET to QCD, lattice 
spacing dependence, light quark mass dependence 
and the quenched approximation. 

The bulk of the matching cancels in the dou- 
ble ratios and the remaining short-distance coef- 
ficients that match lattice gauge theory to QCD 
and HQET to QCD are small and are calculated 
at one-loop level in perturbation theory. Lattice 
spacing dependence is studied using three lattices 
in the range 5.7 < (3 < 6.1. /iai(1) is extrapo- 
lated linearly in to the chiral limit, showing 
a downward trend. The quoted uncertainty is a 
result of including nonlinear terms present in the 
chiral expansion. The chiral extrapolation is in 
fact the largest source of uncertainty. Finally the 
error due to quenching is estimated by summing 
the effect on the running coupling, ag and allow- 
ing an additional 10% uncertainty. 

This is a precision calculation of a phe- 
nomenologically interesting quantity. The re- 
sult agrees well with those from other methods 
(non-relativistic quark models and a zero- 
recoil sum rule |^,|9|). Combining the result 
with measurements by CLEO LEP [|l| and 
Belle H implies 



10^ IK 



cb\ 



45.9 ± 2.4 t}i 
38.7 ± 1.8 
39.3 ± 2.5 ±}:i 



8. rUb 

A new method to determine the b quark 
mass with nonperturbative accuracy was pre- 
sented by Sommer This method has the 
advantage of avoiding the perturbative subtrac- 
tion of power law divergences in the relation 
rhb{^Ji) = ^cont(/^)Afpoic and therefore, taking the 
continuum limit is possible with this method. 
Fig |6| summarises the current status, including 
quenched and unquenched calculations and Som- 
mer's preliminary result. This is in good agree- 



[53]Martinelli-Sachrajda 
[54]GLOK (NRQCD, N,=0) H 
[55]Hombostel et al (NRQCD, N,=0) 
[56]Collins et al (NRQCD, N,=0) 
[57]Heitger-Sommer (N,=0) 
[54]Hombostel et al (N,=2) 
[58]APE (HQET, N,=2) 
[56]Colllns et al (N,=2) 



Average (N,=2) 




4.1 4.2 4.3 4.4 4.5 4.6 4.7 



m,(mj GeV 



Figure 6. Lattice calculations of rah- 



ment with previous determinations. I estimate 
an average of these calculations to be WibiWib) = 
4.30(10) GeV. 

9. SPECTROSCOPY 

In this section I will focus on one quantity: 
the hyperfine splitting in both heavy-light and 
quarkonia systems. Lattice calculations of spec- 
troscopic quantities are generally better con- 
trolled than matrix element calculations of similar 
scope. A major advantage is that no renormalisa- 
tion is required. However, the hyperfine splitting 
(HFS) remains an exception to this rule. 

In the heavy-light sector it is experimentally 
observed that the value of the HFS remains con- 
stant for all flavours. In practice, in the meson 
sector, lattice results are suppressed relative to 
experiment by as much as 20%(40%) for D{B) 
mesons. At this conference Tsutsui described 
a calculation of heavy quark expansion parame- 
ters using NRQCD which includes results for the 
heavy-light meson and baryon splittings. The me- 
son splitting is ~ 30% below the experimental 
value as in other calculations while this discrep- 
ancy does not exist in the baryon splittings. 

It is argued that the discrepancy is a quenching 
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effect since the meson splitting is proportional to 
the wavefunction at the orgin which is suppressed 
in the quenched approximation and to the strong 
coupling which runs differently in the full and 
quenched theories. This however, has not been 
verified by an unquenched calculation of these 
splittings. Lewis Q presented results for the 
charmed baryon spectrum (both QQq and Qqq) 
using anisotropic and NRQCD actions. Despite 
larger statistical errors than in the meson spec- 
trum they see no suppression of the HFS. In fact 
there may be evidence that it is overestimated 
although better data are required to make a con- 
clusive statement. 

Turning to the charmonium and bottomonium 
spectra. Garcia-Perez reported preliminary re- 
sults on behalf of the QCDTARO collabora- 
tion l^ljl . They have calculated the HFS on a very 
fine lattice, /3 = 6.6, using the SW action. They 
take the continuum limit using UKQCD data at 
f3 — 6.0,6.2 and observe significant lattice spac- 
ing dependence. They are currently producing 
their own configurations at a range of f3 values to 
make a more consistent continuum extrapolation. 
New results for the bottomonium spectrum and 
HFS were presented by Manke The calcula- 
tion is done with an anisotropic relativistic heavy 
quark action. Two levels of anisotropy are used 
in the quenched approximation. They also find 
significant scaling violations. 

10. CONCLUSIONS 

Lattice calculations offer the prospect of model- 
independent determinations of hadronic matrix 
elements. The effect of unquenching in calcula- 
tions of leptonic decay constants and their ratios 
is becoming more precise. 

For the semileptonic decay B — > Trliy control 
of the systematic errors in the quenched approxi- 
mation has improved. Problems remain however, 
including simulations with pion momentum above 
IGeV, the chiral extrapolation and unquench- 
ing. An important development is the lattice 
determination of the zero-recoil form factor for 
B — > D*lv. The calculation is systematically im- 
provable and with this method it is conceivable 
that the error on Tb^d*Iu{^) can be lowered to 



under 1%. 
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